Scaling laws and vortex profiles in 2D decaying turbulence 
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We use high resolution numerical simulations over several 
hundred of turnover times to study the influence of small scale 
dissipation onto vortex statistics in 2D decaying turbulence. 
A self-similar scaling regime is detected when the scaling laws 
are expressed in units of mean vorticity and integral scale, as 
predicted in jjj], and it is observed that viscous effects spoil 
this scaling regime. This scaling regime shows some trends 
toward that of the Kirchhoff model, for which a recent theory 
predicts a decay exponent ^ = 10. In terms of scaled vari- 
ables, the vortices have a similar profile close to a Fermi-Dirac 
distribution. 



In recent years, two-dimensional turbulence has re- 
ceived rather large interest because of its applications in 
astrophysics and geophysics and its relative accessibility 
to numerical simulations with respect to fully developed 
three-dimensional turbulence. Two-dimensional flows are 
characterized by the presence of coherent structures (vor- 
tices) which play the role of elementary particles and 
dominate the dynamics. More precisely, the relaxation of 
2D decaying turbulence is a three-stage process: during 
an initial transient period, the fluid self-organizes from 
random fluctuations and a population of coherent vor- 
tices emerges. Then, when two like-sign vortices come 
into contact they merge and form a bigger structure. As 
time goes on, the vortex number decreases and their av- 
erage size increases, in a process reminiscent of a coarsen- 
ing stage. Finally, when only one dipole is left, it decays 
diffusively due to inherent viscosity. 

Two types of studies have been conducted to charac- 
terize this relaxation process: some have focused on the 
precise structure of the vortices (vorticity profiles, u> — ip 
relationships...) while others described how the average 
vortex properties (typical radius, core vorticity, vortex 
number...) evolve with time. There has been some at- 
tempts to predict the final state (the dipole) in terms of 
statistical mechanics of the 2D Euler equation as devel- 
oped in H . It is found that a prediction from the initial 
condition leads to incorrect results due to the effect of 
viscosity which dissipates the high order moments of the 
vorticity during the long evolution of the flow towards 
that state. However, if the constants of the motion are 
evaluated at later times (i.e. before the last merging) the 
prediction gets better and better This implies that 
the statistical theory cannot predict the final state of a 
long viscous evolution but is likely to describe correctly 



the structure of a vortex that forms after a rapid merg- 
ing. It was therefore suggested that the isolated vortices 
of 2D turbulence are sort of quasi-equilibrium states or 
"maximum entropy bubbles" km. In this respect, 2D vor- 
tices share some common features with stellar systems 
like elliptical galaxies which also undergo a mixing pro- 
cess during a phase of "violent relaxation" || . 

Other authors have chosen to disregard the precise 
structure of the vortices in order to study how the typi- 
cal characteristics of the flow evolves with time. In such 
experiments and numerical simulations, it is found that 
the vortex density n, their average radius a and their 
typical core vorticity u> seem to follow power laws. Two 
different scenarii have been proposed. In the first one, 
known as the Batchelor theory M, the assumption of a 
unique invariant (the energy E ~ nw 2 a 4 ), and the oc- 
currence of a unique relevant time scale w _1 ~ t implies 
that the number of vortices n decays as n ~ t~~ with 
£ = 2. This theory also implies the occurrence of a unique 
length scale, the typical distance between vortices nT 1 / 2 , 
of the same order of magnitude as their typical radius a. 
Hence, the total area occupied by the vortices na 2 , or al- 
ternatively the Kurtosis K ~ (ua 2 ) -1 , remains constant, 
implying a ~ t, while the enstrophy Z ~ noj 2 a 2 decays 
like ~ t~ 2 . This is consistent with the selective decay 
hypothesis, which states that for slightly viscous flows, 
the enstrophy is dissipated while the energy remains es- 
sentially conserved. 

However, hyperviscous Navier-Stokes simulations || 
and experiments |J suggest a different scenario in which 
the typical core vorticity lo is an additional invariant. As- 
suming n ~ and the energy conservation such that 
na 4 is now constant, the scaling theory consistent with 
this scenario pj] leads to the slow decrease of the total 
area occupied by the vortices na 2 ~ t~^l 2 (or a ~ t^/ 4 ). 
The enstrophy now decays as Z ~ t~^/ 2 while the Kur- 
tosis increases as K ~ t^/ 2 . The occurrence of an extra 
dimensionless relevant parameter na 2 prevents the de- 
termination of £ from purely dimensional grounds as was 
done within Batchelor theory. 

From the numerical and experimental side, the situa- 
tion is rather confused. Matthaeus et al. (To) performed 
a very long Direct Numerical Simulation of the Navier- 
Stokes equation (DNS) and found that the enstrophy de- 
cays approximatively like £ -1 . More recently, other DNS 
at very large resolution [p"l[ produced a similar decay 
rate Z ~ t~ ' 8 . By contrast, in numerical simulations 
using hyperviscosity (HDNS) |12] the enstrophy decays 



like Z ~ t 
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These hyperviscous simulations show an 
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overall agreement with the second scenario with an ex- 
ponent £ ~ 0.75. 

In a first series of experiments for which 3D effects 
were not fully controlled, Tabeling and collaborators 
jl3| obtained scaling laws compatible with the first sce- 
nario (roughly conserved vortex area coverage), but with 
£ rs 0.44 instead of £ = 2. In a second series of exper- 
iments with stratification the same group obtained 
scaling laws in favor of the second scenario with £ w 0.7. 
In both case, dissipation is provided not via a standard 
viscosity, but mainly via friction at the bottom of the 
experimental apparatus. A simple rescaling however can 
make the experimental system equivalent to a real 2D 
system, with a time dependent viscosity || . 

Theoretical attempts have been made to understand 
and clarify the decay process. Among them, simple mod- 
els describing vortex aggregation process have used point 
vortices following a Kirchhoff-Hamilton dynamics and 
merging via empirical rules derived from imposed con- 
servation laws. The model corresponding to the second 
scenario with constant energy and core vorticity was first 
investigated in leading to £ ~ 0.75 |l^ , |l4| ], in agree- 
ment with HDNS and experiments. Recently, the same 
model was investigated using a renormalization group 
procedure which allows for much larger simulation times 
(3 more decades in time) ||. The true scaling regime is 
only obtained for times much larger than previous simu- 
lation or experimental times, and the asymptotic decay 
exponent is found to be £ ~ 1. Interestingly, in the time 
range comparable to that of HDNS and experiments, the 
function n(t) displays a pseudo-scaling range with an ef- 
fective exponent £ ~ 0.7. An effective three-body theory 
shows that the decay of the total area occupied by vor- 
tices results in a situation where mergings occur princi- 
pally via three-body collisions. A kinetic theory based on 
these three-body processes leads to £ = 1 in agreement 
with the simulations. However, since the conservation 
laws are built in the model a priori, there is a definite 
need for more precise comparisons with DNS. 

Motivated by this observation, we have undertaken nu- 
merical simulations of 2D turbulence at high resolution 
2048 2 , using both normal and hyperviscosity. The goal 
of these simulations was two- fold: first, determine which 
of the two scenarii is more appropriate to describe the 
decline, and whether there is an influence of the numer- 
ical scheme used to dissipate energy; second, determine 
whether there is an asymptotic transition between the 
value £ ~ 0.7 usually reported, and a value £ « 1 pre- 
dicted by the Kirchhoff model, or any other value. This 
imposes to consider a large number of initial vortices, 
so that the decay of their number occurs over several 
decades of time. 

Both viscous and hyperviscous simulations were per- 
formed with a pseudo-spectral code with periodic bound- 
ary conditions. We chose the resolution so that the typ- 
ical size of initial vortices was outside the dissipative 
range. We used a 2048 2 grid for viscous simulations and 



a 1024 2 grid for hyperviscous ones. A random vortic- 
ity field was introduced as initial conditions. The energy 
spectrum corresponding to this initial field is given by 
E{k) = fc 30 /(fc + fco) 60 - Most of the energy is concen- 
trated at the wavenumber fco — 100. This corresponds to 
a situation with approximately 10000 vortices randomly 
localized. The identification of the vortices is based on 
the vortex selection procedure defined by McWilliams || . 
We used the same procedure but with a new definition of 
the vorticity threshold based on both the vorticity maxi- 
mum of the total field and the vorticity maximum of the 
given vortex. 
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FIG. 1. Evolution of the number of vortices and their 
average radius for the three simulations (• and — : DNS; + 
and : HDNS; o and — ■ — : subgrid scale model). 

A good summary of the different scaling laws detected 
in our simulation is provided by Fig. [y. The results of the 
simulation with hyperviscosity mostly confirm the previ- 
ous numerical simulations performed over shorter time 
scales. The number of vortices decays like n ~ t~ - 67 
over two decades in time, while the average vortex radius 
and the Kurtosis increase like a ~ t 015 and K ~ t 30 re- 
spectively. The enstrophy decays steadily over the sim- 
ulation like Z ~ t-° A0 , while the energy remains al- 
most constant, especially at the later stages. Finally, 
the maximum of vorticity is almost conserved, decaying 
approximatively like t~ 0A2 in a first stage, even slower 
(like i" 06 ) in the late stage of the simulation. These 
last scalings are only approximate, since we did observe 
strong local increase of the maximum of vorticity, which 
we associate with strong steepening of the local vorticity 
profile within a few isolated vortices. This is probably 
an artifact of the hyperviscosity, as was previously re- 
ported. Overall, these results are compatible with the 
second scenario with £ w 0.7. 

When normal diffusion is adopted instead of hyperdif- 
fusion, the behavior changes dramatically. Two different 
regimes can be clearly distinguished: in the first one, be- 
tween t = and t w 1 (or equivalently until the total 
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number of vortices has decayed by one order of magni- 
tude), a clean power law n ~ t -0 77 can be observed for 
the total number of vortices, which is close to that ob- 
tained with hyperviscous computations. There is also 
a rather clean scaling law Z ~ for the enstrophy 

within vortices decaying much steeply than in the hyper- 
viscous computations. For any other quantities, a mono- 
tonic but non scaling behavior is observed, with a de- 
crease of the vorticity maximum and of the energy, and 
a slow increase of the average radius and of the Kurto- 
sis. In the second regime (t > 1), rather clean scaling 
laws for most quantities suddenly emerge, and become 
markedly different from the corresponding hyperviscous 
ones. The number of vortices decays like n ~ and 
the average radius increases like a ~ £ - 50 resulting in 
an almost constant vortex area coverage and Kurtosis. 
This regime cannot however be described by the Batche- 
lor theory since £ ~ 1.2 instead of £ = 2 and, in addition, 
u) ~ t~ os and Z ~ i -1 - 3 . 

To test if the observed discrepancy is an effect of the 
finite viscosity, we have also performed inviscid computa- 
tions using a turbulent subgrid scale model described in 
]l5fl . Two different scaling regimes are observed: in the 
early stage of the simulation the Kurtosis remains almost 
constant like in the Batchelor theory but the vortex den- 
sity decays with an exponent £ ~ 0.9 instead of £ = 2. 
At later times, the Kurtosis increases like K ~ t 03 , u> 
becomes nearly constant, and we observe scaling laws 
compatible with the second scenario with £ ~ 0.8. Dur- 
ing all the course of the simulation, the energy is nearly 
constant, like in the hyperviscous case. 

Since all these results were obtained using the same 
initial conditions, they show that the vortex statistics 
is strongly influenced by the dissipative process acting 
at small scale. In particular, the "absolute" scaling law 
exponents, obtained when the quantities are normalized 
with large scale quantities such as the initial energy and 
the size of the box, are not universal. In the spirit of the 
scaling theory, it may however be interesting to deter- 
mine whether a universality class can be detected with 
appropriate JocaJ rescaling of the quantities. The most 
obvious choice is to use as a unit of time the inverse "av- 
erage" vortex vorticity (assumed to be constant in the 
standard scaling theory) defined as Q — \J KZ where K 
and Z are the kurtosis and enstrophy of the vortices. 
As a unit of length, it is then natural to consider the 
integral scale A = u)/>/E built with u> and the energy 
E. These local units are in fact implicitly taken in the 
Kirchhoff model. The scaling laws obtained with these 
units are reported on Fig. | We now obtain a much bet- 
ter agreement between the hyperviscous and the inviscid 
computations, where the scaling laws seem to be com- 
patible with the scaling scenario with £ ~ 0.6, while the 
two regimes of the viscous simulation seem to collapse, 
apart from a small transition zone, into a single regime 
where w/u> ~ (Qt)°, \ 2 n(t) ~ (Qt)' - 6 , a/A - (Cot) - 15 
and Z/{cu) 2 ~ K- 1 ~ {Qt)-°- 33 . 
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FIG. 3. Long-time evolution of the rescaled vortex den- 
sity obtained using inviscid simulations. In the top curve 
the discontinuity arising from a change of subgrid (arrow) 
has been smoothed out. The best fits to the functional form 
n(t) — no/(l+t/to) are displayed. The same data are plotted 
as a function of t+to in the top inset. The effective exponents 
obtained show some trends towards the asymptotic Kirchhoff 
model value £ ~ 1 ^). For short times, coinciding with times 
before the change of subgrid, the data are in agreement with 
the best available Kirchhoff model polydisperse simulations 
(bottom inset). 

The scaling laws obtained in these units are now remi- 
niscent of the early stage of the Kirchhoff simulation. To 
test whether these scaling laws steepen into a regime in 
which £ = 1, one needs to continue the simulations over 
one or two more decades in time, which would represent 
several months of continuous integration using our nu- 
merical resources. This makes longer integrations of the 
viscous or hyperviscous case impossible. However, in the 
inviscid case, we tried to use the flexibility of the subgrid 
scale model to move the large scale/small scale cut-off 
towards larger scales (following the behavior of the in- 
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tegral scale), thereby allowing a gain of computational 
time from 10 to 100. We started the simulation with the 
vorticity field from the DNS at t = 0.3 when the energy 
is small enough at the higher wavenumber. The use of a 
coarser grid induces a small loss of enstrophy and Kurto- 
sis via the filtering procedure used to change the cut-off 
scale in our model. We observed that this change of cut- 
off produces an artificial, cut-off dependent, new scaling 
regime in "absolute" scaling coordinates (using the length 
of the box and the energy as units), but a universal scal- 
ing regime in "local" scaling coordinates (using the inte- 
gral scale and the mean vorticity as units). This universal 
regime is shown in Fig. || for n{t) and clearly suggests 
a local exponent £ effectively increasing with time. The 
extrapolated value for £ is £ ~ 0.79 and £ ~ 0.87 when 
the artificial discontinuity due to the change in resolution 
has been smoothed out (by multiplying the curve after 
the subgrid change by a factor ~ 0.86). 

Our results therefore support the validity of a universal 
self-similar evolution of the vortices for inviscid, or nearly 
inviscid decaying turbulence. This self-similar scenario 
appears universal when appropriate local units are con- 
sidered, and a single exponent in the range £ = 0.8 ~ 0.9 
is found, compatible with that of the Kirchhoff model. 
Finally, viscous effects tend to favor the conservation 
of the vortex coverage area na and modify the scal- 
ing exponents without, however, leading to the Batchelor 
model. 




FIG. 4. Average vortex profile fitted by a Fermi-Dirac dis- 
tribution (full line) and a Gaussian profile (dash-dot). 

We have also found that the vortices present a univer- 
sal profile when the vorticity is normalized by the central 
vorticity and the distance by the typical vortex radius 
defined by the condition w(a) = w(0)/2. This profile is 
represented on Fig. || and has been obtained by averag- 
ing over ~ 30 vortices at different times in the inviscid 
simulation. The error bars indicate to which extent this 
profile can be considered as "universal". As time goes 
on, these bars become smaller showing a trend towards a 
self-similar evolution. We observe that the Fermi-Dirac 



distribution u — <7o/(l + Ae Qr ) provides a very good 
fit to this profile while the Gaussian distribution is less 
accurate (but has only one fitting parameter). In the 
statistical theory of 2D turbulence, the Fermi-Dirac dis- 
tribution maximizes the mixing entropy introduced by 
^ at fixed circulation and angular momentum. Inter- 
estingly, it corresponds to the stationary solution of the 
equation obtained by []l6| from a quasilinear theory of 
the 2D Euler equation. 
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